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Abstract 

I suggest a model based on a softly broken symmetry L e — L M — L T and on Babu's mechanism for 
two-loops radiative generation of the neutrino masses. The model predicts that one of the physical 
neutrinos (^3) is massless and that its component along the v e direction (U^) is zero. Moreover, if 
the soft-breaking term is assumed to be very small, then the vacuum oscillations of v e have almost 
maximal amplitude and solve the solar-neutrino problem. New scalars are predicted in the lOTeV 
energy range, and a breakdown of e-fi-r universality should not be far from existing experimental 
bounds. 



In a model without right-handed (singlet) neutrinos, the three weak-interaction-eigenstate neutrinos 
v e , v^, and v T may acquire |A/| = 1 Majorana masses given by the following term in the Lagrangian: 



— T 



Here, C is the Dirac-Pauli charge- conjugation matrix and M. is a 3 x 3 symmetric mass matrix. One 
may diagonalize M. with help of a unitary matrix U in the following way: 

U T MU — diag (mi, m 2 ,m 3 ) , (2) 

where mi, TO2, and 7713 are real and non- negative. The physical neutrinos i>i, and 1/3 are given by 
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Experiment indicates that two linearly independent squared-mass differences among the three physical 
neutrinos differ by a few orders of magnitude. Indeed, Am 2 tm is of order 10 _3 eV 2 , while Attt-q may be 
either of order 10 _5 eV 2 , in the case of the MSW solution for the solar-neutrino puzzle, or of order 
IO-10 eV 2 , in the case of the vacuum-oscillations ("just so") solution. It is customary to identify ^3 as 
the neutrino which has a mass much different from the masses of the other two, viz., 

12 21 a2„|2 21 _^ I 2 2|_^a2 / r \ 

\ m 2 ~ m l\ = Am < \ m 3 ~ m l\ ^ \ m 3 ~ m 2\ ~ Am atm- ( 5 ) 

Then, the negative result of CHOOZ's search for v e oscillations [Q is interpreted as |f7 e3 | < 0.217, which 
is valid for Am 2 tm > 2 x l(T 3 eV 2 . 
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It has been pointed out S that the assumption of an approximate lepton-number symmetry L = 
L e — — L T (where L e is the electron number, is the muon number, and L T is the tau number) 
may constitute a good starting point for a model of the neutrino mass matrix. Indeed, if there are no 
I Ai| = 2 mass terms then 

rb b \ 

.11= I , (6) 

6 0/ 

where b and r may, without loss of generality, be taken to be real and positive. The mass matrix in 
Eq. (p) yields TO3 = 0, mi = m.2 = b\[\ + r 2 , and 



U = 



r 



V2 
ir 
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v/2(l + r 2 ) v/2 (l + r 2 ) 
1 i 



VT+r 2 
r 



V V 2 (! + ^ 2 ) V 2 ( 1 + 7 ' 2 ) 

This is good for the following reasons: 



vr 



(7) 



1. The negative result of CHOOZ's search for v e oscillations gets explained through U e s = 0. 

2. Since 4 \ U e iU e 2\ 2 = 1, vacuum oscillations of v e with maximal amplitude would occur were mi 7= rri2, 
opening way for the "just so" solution of the solar-neutrino problem to apply. 



It is intuitive to expect r to be close to 1. Now, if r = 1 then v^- 
explains the atmospheric-neutrino anomaly. 



mixing is maximal, and this 



On the other hand, L must be broken, because mi = mi does not allow for oscillations between 
v\ and Vi and a solution of the solar-neutrino puzzle. A good choice, in order to avoid unpleasant 
majorons, would be to have L to be softly broken; this would moreover permit a natural explanation 
for Am,Q <C Am 2 tm . This option has been suggested by Joshipura and Rindani ||; however, in those 
authors' models there is no predictive power for the form of the mixing matrix U, a fact which impairs 
the immediate interest and experimental testability of those models. 

In this paper I put forward a simple model with softly broken L which maintains some predictive 
power. The model is based on Babu's mechanism for two-loops radiative generation of the neutrino 
masses H. I remind that, in general, Babu's mechanism leads to one neutrino remaining massless; 



however, whereas that general mechanism cannot predict the v e 



and v T components of the massless 



neutrino, the specific model that I shall put forward retains the exact-L prediction U e 3 — 0. Moreover, in 
my model there is a rationale for the v e oscillations of maximal amplitude, and for the tiny mass difference 
Am,Q, which allow a "just so" explanation of the solar-neutrino deficit; that rationale is provided by the 
naturalness of the assumption that the term which breaks L softly is very small. 

In my model I just introduce in the scalar sector, above and beyond the usual standard-model doublet 
4> = ( <p + <p° ) , one singly-charged singlet /+ with L — 0, together with two doubly-charged singlets 
g 2+ and h 2+ , and their Hermitian conjugates. The difference between g 2+ and h 2+ lies in that the former 
field has L = whereas h 2+ has L = — 2. The Yukawa couplings of the leptons are L-invariant and are 
given by 



^-{veL <L )\ ji "I ",,L I'L )( )/'/; -I I'rL ~>l I \ ) T H 



m T 
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+/ + [U ("JlC-^l - elC- 1 ^) + f T {v T eL G-\ L - e T L C- l v TL ) 
+e^C- 1 [g 2+ (g^ R + g T T R ) + h 2+ h e e R ] + H.c, 



(8) 



where / T , g^, g T , and h e are complex coupling constants. Notice that, in the first line of Eq. (||), I 
have already taken, without loss of generality, the Yukawa couplings of 4> to be flavor-diagonal; v denotes 
the vacuum expectation value of ip°. 
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The scalar potential V has a trivial part, Vtriviah which is a quadratic polynomial in <jr<j), f f , 
g 2 ~g 2+ , and h 2 ~h 2+ . Besides, V includes two other terms, with complex coefficients A and e: 

v = Vwviai + (\r.rg 2+ + tg 2 ~h 2+ + h.c.) . (9) 



The term with coefficient e breaks L softly. I make the following assumptions: this is the only L-breaking 
term in the theory, and e is small. These assumptions are technically natural in the sense of 't Hooft ||.[[] 

From now on I shall assume, without loss of generality, / M , f T , g T , h e , A, and e to be real and positive. 
Only g^ remains, in general, complex. 

The neutrino mass term M efl does not break L and is generated at two-loops level by the Feynman 
diagram in Figure 1. A similar diagram generates M eT . In both cases, there is in the diagram an inner 
charged lepton which may be either [i or r. It is clear that the mass terms thus generated obey the 
relation 

U /' 

(10) 



M et2 _ ff, 
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Contrary to what happens in Zee's model ||, this ratio of mass terms is not proportional to a ratio of 
squared charged-lepton masses B. As seen before, in order to obtain maximal v^-v^ mixing one would 
like to have r « 1. In the present model, this means that the coupling constants /„ and f T should 
be approximately equal. In Zee's model, on the other hand, one winds up with the rather unrealistic 
constraint / M // T « (m T /m p ) 2 . 

Let us check whether the diagram in Figure 1 is able to yield neutrino masses of the right order of 
magnitude. As we shall see later, we would like to obtain |-M e /J ~ |-A4 e -r| ~ \J Amj tm ~ 10 _2 -10~ 1 eV. 
Now, from the diagram in Figure 1 with an inner r one obtains 

M e ^ = -2\f^f T g T m e m r —-—2, (11) 

(167T 2 ) 



where 



1 ^il dk i2 ,,,2 v > _ ,,,2 I 1 q q 2 - mj q 2 - m 2 . (k - q) 2 - ml 

dy [ y + Xg + l+p , y + x g + x r +p' 
Pin " p In " 

y + x g + 1 - p y + x g + x T - p' 



2(m 2 f -m 2 yo {y + l)(y + x e ) 



+ (1 — x T ) \a.x g + (x T — Xg — y) In x T ] . (13) 
Here, x e = m 2 /m 2 , x T = m 2 /m 2 , x g = m 2 g jm 2 p and 



p = J(y + x g ~l) 2 +4y, (14) 



P = 



y \y + Xg - x T ) 2 + 4yx T - (15) 



The integral in Eq. ( |l3| ) is convergent and may be computed numerically.^] For m e , m T <C rrif and 
Trig Ri m/, one finds / to be of order mj 2 . 

In my estimate of Men 1 shall therefore set / ~ nij 2 . The bounds from e—\x—r universality in [x decay 
and in r decay are f^/rnf < 10~ 4 GeV -1 and f T /mf ^> 10 -4 GeV" 1 |Q; if one allows f^fr/m 2 to be as 
high as 10 -8 GeV -2 , then one obtains 

\M e „\ w 10- 15 A 3t . (16) 

It is reasonable to assume that the Yukawa coupling g T is of the same order of magnitude as the Yukawa 
couplings and / T , and that the dimensionful scalar-potential coupling constant A is of the same order 

1 Notice that the possible L-breaking term / _ / — /i 2 + has dimension higher than the one of g 2 ~ h 2 +, and therefore the 
assumption of its absence is natural. 

2 It is not possible to use the approximations m e = m T = because they lead to infrared divergences. This is not a 
problem, since those divergences are logarithmic and Mefi in Eq. (Jill) also includes a factor m e m T . 
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of magnitude as both mj and m g . This leads to g T /X ~ /u/m/ ~ 10 4 GeV . Fortunately the product 
\g T stays free. In order to obtain |.M eM | ~ 10~ 2 eV it is then sufficient to assume 

A«m 5 «m/ - 10 4 GeV, (17) 
f,~f T ~g T ~ 1. (18) 

Extra factors of order 1 may easily enhance |.M eM | and bring it up to the desired value 0.06 eV. 

The assumption, made in Eq. fll8|), that the Yukawa couplings are of order 1, may seem unrealistic.^] 
However, there are no experimental indications against this possibility when the masses of / + and of g 2+ 
are assumed to be as high as lOTeV.Q For instance, g 2+ mediates the unobserved decay t~ — > pTe + e~; 
however, by comparing that decay with the standard t~ — > [i~i>^v T , one easily reaches the conclusion that 
BR(t~ — > /i~e + e~) should be at least one order of magnitude below the present experimental bound, 
when m g « lOTeV and |g M g T | ~ 1. A more complicated process is e + e~ — » t + t~, which is mediated by 
g 2+ exchange in the t channel. The amplitude A for this process is 

•2 -2 

A = — [v(e)^u(e)] [u(t)7^u(t)] + — — — ^- [^(eh^^e)] [u(r)7 jU 7 5 u(r)] 

" JZ <v(e)r (1 + 7s) «(e)] [u(t) 7/1 (1 + 7s) «(r)] . (19) 



(< — TO 2 



I have used the convenient approximations m e = m T = and sin 2 8 W = 1/4 in writing down the standard- 
model amplitude, and a Fierz transformation in the non-standard contribution. If one defines j = 2m 2 /s, 
z = g 2 / (2e 2 ), and I = 3 (s — to 2 ) /a, then one finds 

da l 2 + l, 9 „v 4 „ Z + l (l + cos6<) 2 , (l + cos6») 2 , . 

oc (1 + cos 2 g) + T cosg + z ■ , , . , J„ a +z 2 ,\ . . —, (20) 



(icos^ I 2 I I \+j + cos6 (l+j + cos( 

where 9 is the angle between the momenta of e~ and of t~ in the center-of-momentum frame. From the 
differential cross section in Eq. ( pp| ) one easily checks that the deviations of both the total cross section 
and the forward-backward asymmetry from their standard-model predictions are completely negligible 
when m g ~ lOTeV, even if g T is as large as 1. 

Except for M. eil and Ai e r, ah other matrix elements of M. break L and, therefore, they will all be 
proportional to the L-breaking parameter e, which is assumed to be small. The matrix elements M.^, 
Mf+r, and M TT arise at two loops from the diagram in Figure 2. In order to obtain a non-zero M ee one 
must go to three loops and use for instance the diagram in Figure 3. In that diagram there are two inner 
charged leptons which may be either /i or t; therefore, there is a contribution to A4 ee proportional to 
m 2 , and that matrix element should not be neglected in spite of it only arising at three-loops level. 

The diagram in Figure 2 clearly leeds to the following relation: 

Mm : Mfi T : Mrr = fl ■ (Ufr) ■ f* = T 2 : r : 1. (21) 
One thus obtains that in the present model 

M = I rb r 2 c rc I , (22) 
\ b rc c J 

where a, b, and c are complex numbers with mass dimension, while r — f^/fr is a real dimensionsless 
number which should in principle be of order 1. The masses a and c are suppressed relative to b by the 
soft-breaking parameter e. 

The mass matrix in Eq. (^2|) immediately leads to two predictions of this model: there is one massless 
neutrino (1^3) and its component along the v e direction, i.e., U e s, vanishes. Indeed, the diagonalizing 

3 Notice however that, in the standard model, the top-quark Yukawa coupling is also very close to 1. 

4 Concerns about the breakdown of perturbativity are only justified for Yukawa couplings >47r, i.e., of order 10 or more. 
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matrix U reads 



U 



costp 
. rsin-0 



—i sin.'i/j 
ia ir cos V> 

VT+7 2 
to * cos ^ 

Vl + r 2 



Vl + r 2 



diag(e 4 <\e lf \l) 



(23) 



cf. Eq. (Q). In the matrix of Eq. (|2^) a = arg [afe* + be* (l + r 2 )] is a physically meaningless phase. 
The Majorana phases 9\ and #2 are necessary in order to obtain real and positive m\ and mi- The sole 
physically observable phase is 2 (6>i — 62) )§)• The mixing angle ip is given by 



tan ip = \J\ + e 2 + 



where 



Ic| 2 fl 



2VT 



(24) 



(25) 



-r 2 |a6* + 6c*(l + r 2 )| 

is a parameter of order e, just as a/b and c/6, and may therefore be assumed to be very small. Thus, ip is 
close to 45°. The amplitude of the vacuum oscillations of v e relevant for the solution of the solar- neutrino 
problem is 4|J7 e if/ e 2| 2 = (l + £ 2 ) , i-e., almost maximal. Thus, the present model favors a "just so" 
solution of the solar-neutrino puzzle. 

The soft-breaking parameter e should be tiny. Indeed, one finds 



Am 2 tm 



\ab* + be* (1 + 7 

\b\ 2 VT+7 2 



(26) 



as we want the "just so" solution for the solar-neutrino puzzle to apply, we must accept e to be of order 
10~ 7 . Such a tiny soft breaking of L may eventually be explained by some new physics at a very high 
energy scale. 

From the non-observation of neutrinoless double beta decay one derives the bound |-M ee | < 0.2 eV 
This is not a problem to the present model. Indeed, as 1TI3 is predicted to vanish, mi and 7712 should 
both be very close to yj Am 2 tm w 0.06 eV. Thus, in the approximation cos 2 %j) = sin 2 tf> — 1/2, one has 



\M e 



(0.03 eV) 



,2 4 (e 1 -e 2 ) 



1 



< 0.2 eV. 



(27) 



Moreover, the phase 2 (9\ — 62) is very close to zero — indeed, it vanishes in the limit of L conservation. 

In conclusion, the model that I have presented in this paper makes the exact predictions 7713 = and 
U e 3 = 0, while it naturally accomodates maximal amplitude v e oscillations and a tiny Am|. Maximal v^- 
v r mixing follows from the reasonable assumption that two Yukawa couplings are almost equal. Neutrino 
masses are small because they are radiatively generated at two-loops level. Indeed, the fact that two 
neutrino masses are as large as 0.06 eV practically forces the new mass scale, at which the extra scalars 
lie, to be in the 10 TeV range; while deviations from e—/j.—T universality in fi decay and in r decay should 
be close at hand. The model requires some physical mechanism for generating a tiny soft breaking of L. 

It is a pleasure to thank Walter Grimus for his generous help and valuable advice, and for reading 
and criticizing two drafts of the manuscript. I also thank Evgeny Akhmedov for calling my attention to 
a potential problem with the model, and for reading the final version. 
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Two- loops Feynman diagram which generates M efl . 
Two- loops Feynman diagram which generates M^ T . 
One of the three- loops Feynman diagrams which generate M. t 
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